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ABSTRACT 

We demonstrate a significant difference in the anguiar momentum transport properties of galactic 
disks between regions in which the interstellar medium is single phase or two phase. Our study is 
motivated by observations of H I in extended galactic disks which indicate velocity dispersions of non- 
thermal origin, suggesting that turbulence in the gas may be contributing significantly to the observed 
dispersion. To address this, we have implemented a shearing-box framework within the FLASH code. 
The new code was used to perform local simulations of galactic disks that incorporate differential 
rotation, self-gravity, vertical stratification, hydrodynamics and cooling. These simulations explore 
plausible mechanisms for driving turbulent motions via the thermal and self-gravitational instabilities 
coupling to differential rotation. Where a two-phase medium develops, gravitational angular momen- 
tum transporting stresses are much greater, creating a possible mechanism for transferring energy 
from galactic rotation to turbulence. In simulations where the disk conditions do not trigger the 
formation of a two-phase medium, it is found that perturbations to the flow damp without leading to 
a sustained mechanism for driving turbulence. The differing angular momentum transport properties 
of the single- and two-phase regimes of the disk suggest that a significant, dynamically motivated 
division can be drawn between the two, even when this division occurs far outside the star formation 
cutoff in a galactic disk. 

Subject headings: hydrodynamics — galaxies: structure — gravitation 



1. INTRODUCTION 

The magnitude of the H I velocity dispersion in the 
outer regions of galactic disks is o bserved to be surpris- 
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ingly consistent, about 7±1 kms ^ (ISantillan et al 
particular exam ples are DDO 154, iDickev et al.l 
ESO_215^?009 ■' Warren 'etalEoM and those discussed 
bv iDib e~a l. 2006 and references therein). Though su- 
pernovae are accepted to be the main energy source for 
interstellar (ISM) t urbulence where star form ation oc- 
curs (see review by lMac Low fc Klessenl[2003 ) the driv- 
ing mechanism for turbulence where this is not the 
cas e is unclear. Theore t ical a rgume nts have been made 
bv iSellwood &: Balbud (|1999l ) and iPiontek k Ostriked 
()2007l ) that this turbulence is m agnetorotationa l insta- 
bility (MRI) turbulence, and by ISantillan et~all ()2007D 
that it is driven through the infall of gas clouds. 

Typically, the stellar disk does not stop immediately at 
a well-defined radius; most commonly the rate at which 
the brightness profile declines switches to a much steeper 
exponential than in the inner disk. This double expo- 
nentia l profile has been characterized by iPohlen et al.l 
()2008l ). Theories predicting a sharp cutof f to star forma- 
tion, such as the Q criteria oflKennicuttI (1989) and the 
cold-phase cri teria of lElmegreen fc Parravano (1994) and 
ISchavd (|2004D do no t naturally explain thi s profil e. How- 
ever, the theory of lElmegreen &: Hunted (|2006f ). which 
includes a second regime of turbulence-triggered star for- 
mation outside the cutoff radius, may p roduce such dou- 
ble exponential profiles. Additionally, iDebattista et al.l 
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()2006( l suggest that an initial single-exponential bright- 
ness profile could evolve into a double-exponential one 
through secular effects. 

Low-level star formation in the extended disk would 
require the cold phase to exist i n at least some regions, 
indeed Ide Blok fc Walter] (|2006[) observed both a warm 
and a cold phase outside the observed limit of star for- 
mation in NGC 6822. 

Self-gravity in conjunction with differential rotation 
has been e xamined as a me c hanism for d riving tur- 
bulence by iWada fc Norlml] (fl999l . [200ll Eoof and 
iWada et all (120021) with global two-dimensional and 
three-dimensional models of a galactic disk including self- 
gravity and cooling. In order to release energy with 
which to drive turbulence, self-gravitational interactions 
must lead to angular momentum transport. Hence it 
is important to analyze the criteria under which self- 
gravitational angular momentum trans por t can occur. 

Following iL ynden-Bell fc K alnaisI ()1972f) and 
ISellwood fc Balbusi (1199911 we introduce the mass- 
weighted average angular momentum transporting 
stress (Txy). The xy subscript results from this being 
evaluated in a local set of Cartesian coordinates around 
one point in the disk. In this frame, x is the radial 
coordinate and y is the azimuthal coordinate. The gas 
perturbed velocity is denoted by u and the gravitational 
velocity is defined as u g = \7^/-/4nUp following 
ISellwood fc Balbui ()1999t ). Then the xy component of 
the stress, which transports angular momentum, is 



{T„j) = {pu 



(1) 
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The first term in this expression is the hydrodynamic 
stress, or TjJ'^'^™, and the second is the gravitational 
stress, or Tfy^'^. The mass- weighted average values of 
these, {Txy), will be used to detect whether a disk has 
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the ability to drive turbulence through the coupling of 
d ifferential rotation to loc al motions. 

IWada fc Normanl (|2007[ ) have performed global simula- 
tions of self-gravitationally driven turbulence in a galac- 
tic disk. In comparison to the initial conditions used 
in this work (see Secti on 13. ip the initial conditions of 
IWada fc Normanl ()2007[ ) lie far into the thermally unsta- 
ble regime. Those simulations are not accompanied by a 
direct analysis of the role of angular momentum trans- 
port in driving the turbulence. This paper carefully ex- 
amines the evolution of a ngular momentum transport. 
iPiontek fc OstrikeJ (|2007t l have performed local simula- 
tions of MRI in the context of a galactic disk with a two- 
phase medium, but do not examine self-gravity. In this 
paper, we present the results of simulations showing the 
difference in self-gravitational angular momentum trans- 
port properties of single-phase and two-phase galactic- 
type disks. We also demonstrate the onset of this mode 
of angular momentum transport when a single-phase ac- 
creting disk gains sufRcicnt surface density to form a two- 
phase medium. This was done by a series of models of 
local disk patches, with either fixed mass or continuous 
accretion. The possible regimes of behavior for a local 
patch of galactic disk are then discussed. 

2. METHODS 

We began with the standard FLASH 2.3rl code 
(jFrvxell et all |2000D which, among other features, in- 
cluded support for parallel adaptive mesh hydrodynam- 
ics and self-gravity in periodic boundary conditions. We 
implemented support for shearing-box boundary con- 
ditions and fictitious forces, shearing-box self-gravity, 
and modified the hydrodynamics for the shearing frame. 
The background shearing flow was split from the hy- 
drodynamics, in a manner similar to 'Masset' (2000) and 
[Gammi c (2001i). The e volution of linear she aring waves 
has been t e sted, as in iJohnson fc Gammi^ (j2005l ) and 
iShen et all (l200l . 

To handle the combination of both self-gravity (Pois- 
son's equation) and the shearing frame (Hill's equations; 
iHilll . |l878) three choices are made. First, the vertical 
background force term, originally fl'^z, was transitioned 
to in the range Zc < z < 2zc, where Zc is a param- 
eter, by means of a cubic interpolation. The vertical 
external gravity was hence zero at the vertical bound- 
aries. Second, the vertical fluid boundary is set to either 
zero-gradient (outflow) or a fixed- value (inflow). Third, 
self-gravity is solved with a periodic boundary condition 
in the vertical direction, and the simulation box kept 
tall and thin so that the combination of the self-gravity 
and the vertical background force produces a boundary 
condition essentially equivalent to an isolated boundary 
condition for the disk. 

The shearing box is specified by two parameters, the 
angular frequency CI and the radial rate of change of this 
quantity, q = —dlnCl/dlnR. To simulate an outer galac- 
tic disk, these are chosen as £7 = 8.9 x 10^^® s^^ and 
q = l. 

As the shearing box is a local model, global self- 
gravitational phenomena such as spiral arms and bars 
cannot be captured. The shearing box also forces the an- 
gular momentum trans port to be purely local, i.e . within 
the scale of the box (Balbus fc Papaloizoul[l999l ). In the 
context of this work, this is an advantage, as the mecha- 



nisms produced will be local and minimal in nature, not 
relying on the large-scale structure of the galaxy, and 
hence are, perhaps, more generally applicable. 

2.1. Cooling 

We investigate the significant difference in behavior 
between single- and two-phase disks. The formation of 
a two-phase disk is dependent on the radiative cooling 
used. To include the local effects of various optically 
thin limit cooling and heating mechanisms a net heat 
loss function is defined as £(T) = pA(T) — F, where 
F is a constant following [Brandenburg. Korpi. fc Med 
(|2007| ) and others. We use a heating par ameter F = 
0.015 ergg~^ s~^ as in [Brandenburg et al.l p007.) and 
IPiontek fc Ostrikeil (l2007l). The form of A(r) is a piece- 
wise power-law fit from " iBrandenburg et al.l (|2007t ). The 
resulting heating equilibrium curve is shown in Figure [TJ 
For calculating temperature, the mean particle mass is 
fi — 0.62 and a 7 = 5/3 ideal gas equation of state 
is used. To ensure that the Jeans l ength is resolved, 
the Je ans length limiter is taken from IWada fc Normanl 
(|2007( ). This is designed to prevent artificial fragmenta- 
tion by enforcing a minimum value of the Jeans length 
(|Truelove et al.l[l997i) . The Jeans limiter artificially halts 
collapse where is can still be resolved without artificial 
fragmentation. Runs were stopped when the computa- 
tional expense became prohibitive. Whether the collapse 
stops in real disks, and how, cannot be addressed by these 
models. 

3. RESULTS 

3.1. Fixed Mass Disk 

The important parameter in fixed-mass disk models 
is where the phase of the mid-plane lies in relation to 
the critical curve for thermal instability. Disks in which 
the densest parts have perturbations that are sufficiently 
strong to trigger thermal instability, form a two-phase 
medium. Seven models with disks of fixed mass and cool- 
ing were run. The models had parameters such that the 
mid-plane conditions varied around the heating-cooling 
equilibrium curve on the p-P phase diagram. These mod- 
els provided characterizations of the possible two-phase 
and single-phase disks. The parameters of these runs 
were a box 1.2 kpc x 1.2 kpc x 4.8 kpc, with no minimum 
level of refinement. The standard resolution run used 
four levels of refinement at the mid-plane or a mid-plane 
resolution of 128^ and 9.375 pc. The initial condition was 
a hydrostatic equilibrium isothermal disk, with 1% ran- 
dom density perturbations added to each cell. Runs with 
surface density E = 6 Mq pc~^ and initial temperature 
2500 K, the standard mid-plane resolution, and double 
the standard resolution developed a two-phase medium. 
An additional variation with S — 15Mqpc^^ and ini- 
tial temperature 10, 000 K also developed this two-phase 
medium. The unstable filaments connecting the lumps 
are gradually accreted into clouds and an unstable layer 
exists as a sheath around the cool clouds. A Q-stable 
disk of thermally stable gas remains extending far above 
the clouds. The phase diagram of this two-phase medium 
consists of points lying on the warm thermal equilibrium 
line, a scatter of points proceeding across to the cold equi- 
librium line below the peak warm equilibrium pressure, 
points lying on the cold equilibrium line, and a small 
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Fig. 1. — Run Al gas phase distribution at t = 2.73 X 10^® s with 
equiUbrium curve and temperature limits 
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Fig. 2. — Run Al stress - showing the transition from stable 
to unstable evolution as the cool phase develops, the gray vertical 
stripe marks the onset of the formation of the cool phase, the jump 
in (Tjy^^) is approximately a factor of 10* 

number of points lying on the Jeans limiter line. This is 
similar to the phase diagram for the late time of the ac- 
creting disk shown in Figure [TJ When the disk is single- 
phase, it evolves toward a smooth shearing flow with- 
out significant angular momentum transporting stress. 
The four remaining runs evolved to this state. Addition- 
ally, Toomre-stable isothermal gas disks (without cool- 
ing) also behaved in this manner. 

Even in the two-phase disk runs, the tot al mass- 
weighted valu e of the Toom re Q parameter (Toomr^ 
[1964; Goldrei ch fc Lvnden-B ell 1965) stays weU in the 
range of stability. For example, Q > 4.5 for all times in 
the first two-phase run described, even though the cool 
phase by itself is Toomre-unstable. The self-gravitational 
stresses of two-phase fixed mass disks are essentially the 
same as the two-phase part of the accreting disk to be 
discussed in Section [221 

3.2. Accreting Disk 

Run Al starts with the grid as in the standard res- 
olution fixed-mass disk in the previous section, but 
with an initial condition of uniform very low den- 



sity (10~^* gcm"-^) gas. Inflow z-direction bound- 
aries were specified with the inflow rate set as 1.146 x 
10~® M0 pc~^ yr. Along the 16 x 16 grid of cells at the 
z-boundary, the density is 1% greater on odd- numbered 
rows. This provides a minimal but nonangular momen- 
tum transporting perturbation to the flow, so that the 
imposed perturbations will be larger than those aris- 
ing from purely numerical effects. With these boundary 
conditions, the disk built up and eventually reached the 
critical temperature/pressure for thermal instability and 
formed a cool phase. The inflow rate was chosen to be 
slow so that the cool phase formation is not triggered 
directly by the inflow pressure. 

The disk starts at very low surface density and slowly 
builds. As the infalling gas from either boundary collides 
with the disk it drives turbulence in the disk as it builds. 
This is however purely a function of the inflow boundary 
conditions. As the disk surface density grows, so does 
the mid-plane pressure, until eventually it is sufficient 
to trigger thermal instability, and the cool phase read- 
ily forms. Turbulence driven by the infalling gas creates 
density perturbations which the formation of the cool 
phase amplifies. As these cool lumps form and rapidly 
accrete gas, the gravitational interactions between them 
lead to vastly increased angular momentum transport- 
ing gravitational stress, as shown in Figure [21 When 
the gravitational stress amplifies as a result of the phase 
transition the change is dramatic: in the case of run Al 
it is 4 orders of magnitude. After the cool phase forms, 
the phase diagram is like that displayed in Figure [l] Gas 
is introduced at a point near the heating-cooling equilib- 
rium, initially expands and cools as it falls in, then re- 
compresses and heats, rising above the equilibrium line 
toward higher densities. There is then a layer of gas lying 
along the equilibrium line, until the critical pressure for 
the phase change is reached. A further scatter of points 
on the phase plot shows the layer of cells around the cool 
clouds where gas is unstable and cooling onto the clouds. 
Finally, points denoting the cool clouds are distributed 
along the cool equilibrium line. 

The run A2 repeats the same setup as run Al, but 
with half the mass inflow rate. A2 has the same pattern 
of rapidly growing stresses, with the transition at ap- 
proximately the same surface density, demonstrating the 
critical importance of surface density rather than sensi- 
tivity to the imposed inflow rate. Between runs A2 and 
Al the pre-transition hydrodynamic stresses vary in a 
different manner, demonstrating that they are indeed a 
function of the imposed boundary conditions. 

4. DISCUSSION AND CONCLUSIONS 

Two primary conclusions can be drawn from this work. 
First, thermal instability can trigger self-gravitational 
angular momentum transport. Second, the development 
of the cool phase leads to rapidly growing local self- 
gravitational angular momentum transporting stress. 

The thermally stable cooling disks of fixed mass de- 
cayed toward a smooth flow state. Similarly, before the 
onset of thermal instability, the accreting disk models did 
not have significant self-gravitational stress. However, 
the combined results of the fixed-mass disk with cool- 
ing and the accreting disk simulations show that ther- 
mal instability and the formation of a two-phase medium 
can provide a mechanism for producing enhanced angu- 
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Fig. 3. — Three regimes of a galactic disk, the P — p insets 
show which parts of the phase diagram are populated around the 
cooling-heating equilibrium curve 

lar momentum transport by self-gravitational stresses. 
This will have two consequences: the galaxy's rotational 
energy is converted into local motions giving rise to ve- 
locity dispersion in the gas; and it enhances the rate at 
which mass, primarily in the form of cool gas clouds, 
is transported inward enhancing the surface density of 
more central regions of the galactic disk. Increased sur- 
face density may in turn lead to star formation. 

If the observed roughly constant velocity dispersion in 
the outer galactic disk is to be driven by a local self- 
gravitational mechanism, the simulations in this work 
suggest that the cool phase must exist everywhere in the 
disk where the velocity dispersion is seen. However, as 
the surface density of the disk drops, it will be increas- 
ingly difficult to achieve the critical pressure for the ex- 
istence of the cool phase at the disk m id-plane. Detailed 
observ ations, extending the work of Ide Blok fc Walter! 
(|2006f ). might settle this by directly detecting a cool 
phase in extended, low-surface-density galactic disks. 

Finally, a picture of the phase structure and dynamic 
structure of the outer regions of a galactic disk can be 
drawn. Three regimes, as shown in Figure [3l can be 
identified, with the divisions between them drawn by the 
cutoff in star formation and the gas single-phase to two- 
phase transition. 

I A region of single phase gas, closely following the 
heating-cooling equilibrium curve on the phase di- 
agram, with Q > 1, and with pressure too low to 
support a two-phase medium. Includes small per- 
turbations. The total disk is Toomre-stable. 

II A region where pressure is sufficient for cool clouds 
to form from existing small perturbations, forming 
a two-phase medium, enhanced angular momen- 
tum transporting stresses result. The total disk 
is Toomre-stable. 

Ill The stellar disk, where the conditions for star for- 
mation are widely satisfied, and feedback (super- 
novae and outflows) drives turb ulence. The tota l 
disk has Toomre Q satisfying a iKennicuttI (|1989l ) 
relation. 

Regions I and II are simple to form - thick warm gas 
disks are physically simple, and the stellar disk is well 
known to exist. There are already a number of elements 
pointing to the nature of region II. iPiontek fc Ostriked 



(|2007D point out that the lSchavd (|2004D hypothesis, that 
the formation of the cool phase necessarily implies star 
formation, fails based on current observations showing 
that the phase change can occur outside the cutoff of 
star formation. That a region where turbulence must 
be driven by nonstellar proces ses while the co ol phase 
exists has also been argued bv lOstriked (|2007[ ). In the 
same work, they also find that magnetohydrodynamic 
effects in a region with cool clouds, such as in region 
II, may suppress star formation. The presence of star 
formation sets the division between regimes II and III. 
This cutoff may not be sharp and it could be expected 
that locally star formation may happen in the cool phase 
of region II. The presence o f stars, and their feedback 
(particularly supernovae, see i Mac Low fc Klessenl . [200l 
changes the dynamics and phase structure of gas from 
II to III. Between regimes I and II the dividing line is 
drawn on the existence of the cool phase. In terms of 
the importance of local self-gravity in galactic disks, this 
work has shown that there is a significant energetic divide 
between regimes I and II. 

The self-gravity and hydrodynamic stresses presented 
here in media with nonzero shear can be calculated and 
can be of similar magnitudes in nonshearing cases. A ver- 
sion of run Al without shear was preformed to demon- 
strate this. The definition of the stresses does not de- 
pend explicitly on the existence of nonzero shear. How- 
ever, without shear, these stresses are not significant as 
they cannot extract energy from the background shear or 
transport angular momentum without a nonzero shear to 
couple to. 

In regime II simulations in this work, the mass 
weighted average power per unit volume from self- 
gravitational stresses is typically on the order of 
10"'^^ erg c m~^ s~^, a s ignificant energy input in such 
a context ()Mac Low fc Klcssen 2004). We have also 
demonstrated that the presence of self-gravity leads to 
larger velocity fluctuations as follows. In a fixed-mass 
disk, we evolved a S = 6 pc^^ disk with an initial 
temperature of 1250 K until a two-phase medium had 
developed. At this time the mass wighted RMS velocity 
of gas below 1000 K was 1.3 km s^^. This state was then 
used as input for two further simulations: one with the 
normal self-gravity, and the second with a fixed vertical 
gravity that was obtained by horizontally averaging the 
vertical self-gravity forces in the input state. When it 
became computationally infeasible to evolve the simula- 
tions further, the self-gravity case had a mass wighted 
RMS velocity of gas below 1000 K that had grown to 
2.6 km s~^ and was climbing, while the same quantity in 
the fixed-gravity case had fallen to 0.8 km s^^ and was 
declining. 

The extent to which extended H I disks have regions of 
type II must be determined by observations detecting or 
showing the absence of a cool component. One existing 
example is the observations bv lde Blok fc Walter] ()2006f ) 
of NGC 6822, which suggest that part of this disk may be 
characterized as regime II. Even in the relatively small 
box used in Section [XTl the mass wighted RMS velocity 
of gas below 1000 K is 2km s~^. This is in the same 
ballpark as the observed H I velocity dispersions. 
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